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ABSTRACT 

We begin a program of work aimed at examining the interior of a rotating black 
hole with a non-zero cosmological constant. The generalisation of Teukolsky's equation 
for the radial mode functions is presented. It is shown that the energy fluxes of scalar, 
electromagnetic and gravity waves are regular at the Cauchy horizon whenever the surface 
gravity there is less than the surface gravity at the cosmological horizon. This condition is 
narrowly allowed, even when the cosmological constant is very small, thus permitting an 
observer to pass through the hole, viewing the naked singularity along the way. 
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1. Introduction 

The fate of an observer who falls into a black hole is an interesting and important issue 
within the context of general relativity. The situation of interest is when the black hole is 
charged or rotating, because then the spacetime can be continued beyond the black hole 
into another universe, with naked singularities and causality violation to contend with. 

We know that for a black hole in empty space the observer's progress is impeded 
by a singularity at the Cauchy horizon [1,2], but the situation is different if there is a 
cosmological constant. In the non-rotating case we now know that there is a finite range 
of parameters for which the Cauchy horizon is fully stable [3-6] . In these cases an observer 
could pass into the hole and see the naked singularity. Since a cosmological constant does 
not contradict the laws of physics in any known way, this journey constitutes at the very 
least a valid thought experiment in which the predictive power of general relativity breaks 
down. The evolution of the spacetime beyond the Cauchy horizon requires a knowledge of 
'initial' conditions at the singularity and this lies beyond the limitations of current theory. 

In this paper we examine the stability of the Cauchy horizon for a rotating hole 
when a cosmological constant is present by consideration of the energy fluxes of scalar, 
electromagnetic and gravity waves propagating on the black hole background. We begin 
by separating the field equations into equations for radial and angular functions. For 
this step in the calculation the introduction of a cosmological constant becomes an asset 
because it makes manifest the symmetry of the metric under the interchange of the radial 
coordinate r and the coordinate fx = a cos 9. It is the strength of this symmetry combined 
with the algebraic character of the curvature tensor that makes possible the separation of 
the fields and ultimately makes the problem solvable. 

The significance of gravitational wave equations to the stability lies in the fact that 
they result from perturbing the full set of Einstein's field equations [1]. Therefore an 
arbitrarily small perturbation, provided it does not grow in size anywhere in the spacetime 
region of interest, can be analysed by purely linear equations. This is to be contrasted 
with the unstable case, familiar from the ordinary Kerr metric, where it is necessary to 
solve non-linear equations because an initially small perturbation becomes large near to 
the Cauchy horizon [2]. 

As a matter of fact, gravitational wave propagation is only part of the full gravitation 
perturbation analysis because there are other degrees of freedom in the metric, but it 
ought to contain all of the physical degrees of freedom relevant to the stability question. 
To be certain, it only remains to be shown that the remaining gravitational perturbation 
equations separate into similar wave equations to deduce stability in the full gravity-matter 
system. We hope to establish the complete separation of the equations in a later paper. 

We find stability requires that the surface gravity k of the cosmological horizon should 
be larger than the surface gravity of the Cauchy horizon; the same condition as in the 
charged case. The parameter range for which this condition holds are found in the following 
section, in which various properties of the spacetime are explained. The parameter range 
is small, but finite. 

2. The spacetime 

The metric of a stationary black hole in de Sitter space was first given by Carter [7] 
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and takes the form, 



ds 2 = p 2 (A r 1 dr 2 + A ^dp 2 ) + p 2 A fl u) 1 ®u) 1 -p 2 A r u 2 ®u 2 , 



(2.1) 



with one-forms, 



X 00 = dt — a~ a r d(j), \ u = dt — a~ a d(p. 



(2.2) 



The metric is parameterised as follows, 



(a 2 + r 2 )V2, 
(a 2 -^ 2 ) 1/2 , 



A r = ( a 2 + r 2 )(l - ±Ar 2 ) - 2Mr + Q 2 
A M = (a 2 -^)(l + |A/, 2 ) 



(2.3) 



and 



p 2 (r) = r 2 + p 2 , 



X 2 = l + ^Aa 2 , 0(r)=a/(r 2 + a 2 ). 



(2.4) 



The function O is explained below. Note that the subscripts are used to denote functions 
that depend only on the coordinate indicated by the subscript, and a prime will always 
denote differentiation with respect to this coordinate. 

The metric is a solution of the Einstein-Maxwell equations with a non-zero cosmo- 
logical constant A and the electromagnetic vector potential of a charge Q, 



To see that this metric really does represent a black hole we introduce the angle 9, where 
p = a cos 9. The surfaces of constant radius r and time t are distorted spheres, parame- 
terised by spherical polar coordinates 9 and 4> in the usual way. 

The metric has coordinate singularities at the roots of A r = 0. For a range of pa- 
rameters shown in figure 1 there are four real roots, which shall be labelled in decreasing 
order by r±, r2, rs and ta. The roots T2 and rs represent the outer and inner horizons 
respectively of the black hole, whilst r\ represents a de Sitter horizon. The fourth root ta 
is negative and non-physical. 

It is possible to set up coordinate systems for which the metric is regular on any 
particular horizon. Following section 56 in ref. [1] , in the region outside the event horizon 
there exists two simple families of null geodesies with tangent vectors 



where a dot denotes a derivative with respect to the geodesic affine parameter and u> is 
a constant. Corresponding 'light-like' coordinates which are fixed along each family can 
therefore be defined by u = t — r* and v = t + r* , where 




(2.5) 




(2.6) 



r 



dr* = X 2 K 1 °rdr. 



(2.7) 
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Additionally, for rotating black holes, there is a unique combination of the killing 
vectors that is null on the (future) event horizon, / = dt + 0(^)^0. This leads us to 
introduce an azimuthal angle coordinate 

2 = 0-Q(r 2 )t (2.8) 

on the horizon, which is constant along the integral curves of /. Physically, this corresponds 
to the dragging of inertial frames near to the hole. 

The null-killing vector on the horizon can be used also to define the surface gravity, 
V I 2 = Kil. The explicit expression for Ki follows from squaring both sides of the expression 
and taking the limit as r — > ri [8] , 

Ki = ^A X - 2 a- 2 l[\r i -r j \. (2.9) 

(There is a typographical mistake in ref. [9]). 

The extension of the metric through the event horizon was given by Carter in ref [7] 
(with some corrections given in ref. [8]). The regular coordinate system is given in terms 
of the surface gravity, 

U = -K 2 - 1 e" K2M , V = K^e K2V (2.10) 

The coordinates u and v will be defined inside the event horizon as shown in figure 2. 

It is possible to extend the metric through each of the horizons in turn and obtain the 
complete Penrose diagram, a part of which is shown in figure 3 [7,8]. The fully extended 
spacetime stretches arbitrarily far in each direction. It is possible to find world-lines that 
pass through the black hole from region /, through regions II and III and on to the 
outside. 

The exterior geometry of a collapsing star would only take in part of the picture, the 
remainder being replaced by the geometry of the stellar interior. In order to retain the 
axial symmetry of the solution it is also necessary to place a second body of equal mass 
around the point of the universe that is antipodal to the centre of the star. This second 
body can be thought of as representing the rest of the matter in the universe. 

Inside region III lies the naked singularity and the possibility of causality violation 
from closed timelike curves. These curves pass beyond the singularity at r = 0, which only 
appears in the // = plane, into the next region in the Penrose diagram. The problem for 
an observer who tries to enter these regions lies in crossing the Cauchy horizon between 
regions II and III at r = r 3 . 

As an observer crosses the Cauchy horizon any radiation from outside appears to be 
either red-shifted or blue-shifted depending on the relative values of the surface gravity at 
the Cauchy and de Sitter horizons, red-shifted if the Cauchy horizon value is the smaller. 
In the non-rotating case, this was also the Cauchy horizon stability condition [4,5]. 

It is possible to split up the parameter space of the black hole metrics into regions 
depending on the relative sizes of the surface gravity on various horizons. We write the 
difference between the surface gravities at ri and rj as 
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where pij is a polynomial in the radii. 

If two of the radii are equal, then the surface gravities of the respective horizons both 
vanish. In other cases the relative magnitudes of the surface gravities are determined by 
the signs of the polynomials. Explicit expressions for these polynomials are, 

P12 = (rir 2 + r 3 r 4 - 2a 2 ) (r 2 - r|) 

P13 = (*"ir 3 + r 2 r 4 )(r 1 + r 3 ) 2 + 2(nr 3 - a 2 )(rir 3 - r 2 r 4 ). 

The black-hole parameters for which p^ = can be found by first writing p^j in 
terms of x = r\r 2 — r^r^ or x = rir 3 — r 2 r 4 in the case of pi 2 or pi 3 respectively. From the 
conditions on roots of A r = it is possible to deduce quite independently that x satisfies 

M 2 = ±A 2 x 3 + ^A(l - iAa 2 )x 2 + ±A(a 2 + Q 2 )x + (1 - ±Aa 2 )(a 2 + Q 2 ). (2.13) 

Solving pij = for x gives an expression x(a,Q,A) which can be substituted in this 
equation to give a condition M(a,Q, A) for equal surface gravities. 

The situation in which the surface gravities on the de Sitter and event horizons are 
equal gives rise to gravitational instantons and was discussed in reference [8]. From equa- 
tion (2.12), the solution for this case is x = 2a 2 and leads to a condition from (2.13), 

M 2 = Q 2 X 2 + a 2 x 4 . (2.14) 

The borderline of the stability criterion is that K\ = re 3 . This can be found by solving 
Pi3 = numerically and substituting the result into (2.13), leading to the condition shown 
in figure 1. The line lies very close to the line where the inner and outer horizons coincide. 
The region of stability between the two lines OA is small, but finite in a uniform measure 
on the space of parameters M, Q and a for fixed A. 

3. The scalar wave equation 

We will take a massless scalar field with wave equation 

V 2 $ = 0, (3.1) 

to demonstrate how the wave-equation separates on the black hole background [10]. 
In terms of coordinates the Laplacian becomes 

v 2 = ^a r A r a r + + ^d 2 - ^-d 2 (3.2) 

P 1 p 2 A M A r 

where d\ and <9 2 are dual vectors to ui\ and u>2, 

2 2 

di = --^{<y\ dt + adcf,) d 2 = -^(a 2 d t + ad^). (3.3) 

We look for separable solutions of form 

$ = R(r)S(ij)e- lujt e tm *. (3.4) 
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The angular and radial terms in the wave equation separate immediately, with angular 
part 



We require that the solutions be regular at /j, = ±a. Scaling /i by a shows that 
the eigenvalues A depend upon the parameters m, ua and Aa 2 . We also add another 
parameter / to label the position of the eigenvalue in sequence, starting from the value \m\. 
The corresponding eigenfunctions will then be denoted by S(lmuj; y). Some values of the 
angular eigenvalues obtained numerically are tabulated in tables 1-3 for various parameter 
values. The cosmo logical constant has been re-expressed there as A = 3/a 2 . 

For A = the equation reduces to the equation for spheroidal wave functions, 



The spheroidal eigenvalues denoted by \™ [11] are then related to A by a shift of 2amu. 
We shall be interested later in the analyticity properties of the angular functions in a 
strip of the complex uj plane below the real axis. Numerical studies for spheroidal wave 
functions indicate that there are branch cuts away from the real axis and suggest that 
these functions are regular in the region \uj\ < cl 2 for some constant c [11]. 

For small values of A it is possible to relate the properties of the angular functions to 
the spheroidal wave functions by a suitable approximation scheme. The angular functions 
are the stationary points of an integral, 




(3.5) 



where 




(3.6) 



((a 2 - y 2 ) S')' + [A + 2amu - u 2 (a 2 - y 2 ) - a 2 m 2 (a 2 - y 2 )' 1 ] S = 0. (3.7) 




(3.8) 



We can separate the A = part and write 




(3.9) 



n 



where 




[S U S 2 ] = P dy {y 2 alS[S 2 + (1 + a 2 /al)( 



am — a^uj) 2 S1S2] 



(3.10) 



J — a 



am — a 2 oj) 2 SiS2- 



These give series for small A, 




(3.11) 



n 



n 
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where is the spheroidal wave function. Thus 

A (D = j(D [5 (0) j S {0) ]t 5 (D = £ (A (o) _ A (0))-ij(D [5 (0) j 5 (0)] 5 ( 0)- (3 . 12) 

This result was used to check the values in the tables. Continuing further with the series, 
it is possible to show that the are regular functions of u when is regular. 
Returning now to the wave equation and taking the radial part, 

K 2 

-d r (A r d r R) + XR - —^R = 0, (3.13) 

where 

^ r = X 2 (am-(7 r 2 w). (3.14) 

It is possible to write 

a r 9 r (A r 9 rJ R) = /* 2 ) (^) (3-15) 

This allows the radial equation to be written in the form of a scattering problem, 

^2 + V(r)j OrR = (3.16) 
where the potential V(r) takes the form 

V = cu — mil) -. — 7 -. — 5- — — . 3.17 

This potential has been plotted as a function of r* between the outer and inner horizons in 
figure 4. There is a potential well in the region covered by the plot and a potential barrier 
outside of the event horizon. The main difference from the non-rotating case is that the 
values of the potential as r* — > oo and r* — > — oo need no longer be equal. 



4. The electromagnetic wave equations 

In this section we will examine the Maxwell field equations on a Kerr-de Sitter back- 
ground using the Newman-Penrose [12] formalism. This development follows the lines of 
the analysis for a metric with vanishing cosmological constant [13-15]. The main difference 
lies in exploiting to the full the symmetrical appearance of the metric with respect to r 
and [i. 

We begin with an orthonormal tetrad of vector fields, 



e„ = 



r -d r , e 2 = — 7 ={a r d t +ad (j} ) 



-d^, ei = 




(4.1) 



(a d t + ad < j t ). 
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The vectors dual to uj 1 and uj 2 have been used as in section 2. Newman-Penrose formalism 
consists of writing the connection forms in a complex null-tetrad basis, 

I = -7^(e r + e 2 ), n = -4=(e r - e 2 ) 

i i ( 4 - 2 ) 

m = ^(e M + zei), m = ^(e^-«ei). 

These vectors satisfy /.n = —1 and m.rh = 1. (Inside the event horizon we will use the 
modulus of A r for scaling the tetrad vectors and choose / always to be in the direction of 
dy, where V is the Kruskal coordinate.) 

The notation for the covariant derivatives along the tetrad vectors consists of 

D = V h A = V n , S = V m , 5 = Vrn. (4.3) 

The connection forms are then expressed as in the following table, 



Dl = (e + e)l — Rm — Km, Dm = (e — e)m — ku + nl, 

Al = (7 + 7)/ — fm — rm, Am = (7 — 7)777. — rn + Pi, 

dl = ((3 + a)l — pm — am, dm = ((3 — a)m — an + XI, 

61 = (a + (3)1 — am — pm, dm = (a — (3)m — pn + pi. 



(4.4) 



For the Kerr-de Sitter metric some of the connection components vanish, 

K = a = X = u = 0. (4.5) 

The vanishing of these connection components is a feature common to all black hole solu- 
tions. The remaining connection components are 

1 c r— 1 c r— 1 a/AT' 

P = P = E-TV A r> e = 7 = -r^^V A r + 



V2p 3 v r ' ' 2v/2P 3 V 2v/2 P ' 

* C /Hi— A, 7 C r-r— 1 a/AT 



(4.6) 



v/2> 3V M ' 2v/2P 3V " ' 2V2 P 

where 

c = r + ip. (4-7) 

The electromagnetic field is described by three complex scalars <po, (pi and <p 2 , with 
the radial fields in 0i and the axial fields in 0o + 02- They satisfy Maxwell's equations, 

D0i - 50o = (tt - 2a)4> + 2p0i - K<j) 2 , 

5(pi - A0 O = (p - 27)^0 + 2r0i - a<j) 2 , ^ 
D(f> 2 - 5(f>i = (p - 2e)0 2 + 2n(j) 1 - A0 O , 
5(p 2 - A0i = (r - 2/3)^2 + 2/70i - W>o- 
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We follow reference [1] and introduce new radial and angular derivatives defined by 



V r 



%k t a; , iK r a; 
+ aT + n A? ^ = ^" aT'+^aT' 



and 



K 



A' 



A', 



- n "^ + A„. +n A' Ln ~^ A,,. A„. ' 



(4.9) 



(4.10) 



with K r and denned in equations (3.14) and (3.6). These are respectively purely radial 
or angular in character. Useful properties of these derivatives are that 



n 



V n + - = (cT^ n (cp, V n = A~ n V A? 



(4.11) 



and similarly for C. For scalar functions that have the dependence expi(m</> — ut) on <p 
and t it is possible to write the directional derivatives along the tetrad vectors as 



1/2 



D= S 2>o, A = - ^ Pj, 5 = 



2p 2 



1/2 



2p 2 



1/2 



(4.12) 



We can now write Maxwell's equations in terms of V and £, using the connection 
components from equation (4.6). We will rescale the Maxwell fields first, to ipi = c(fii, and 
then the equations become 

^A r (v + ^ (f! =^/A^ (c\ /2 + £j if 
1 



'A r V 



'1/2 " 1 \ ¥0 



= a/AT I Cq - - )(p 1 



A r V 1/2 



<P2 



\ [ 4 - j 



(4.13) 



The field </?i may be eliminated from (4.13), 



A„ Co 



>t 1 



A^ £ T 1/2 + z }<fo = -VA r £> + 



1/2 



V? (4.14) 



From the forms of T> and £, and of _KT r and K^, this can be simplified to 
(£_ 1/2 A/! /2 + P_ 1/2 A r P} /2 + 2zw X 2 c) </?o = 0. 

Similarly, 

(£i 1/2 A M £ 1/2 + Pl 1/2 A r P 1/2 - 2zw X 2 c) </? 2 = 0. 



(4.15) 
(4.16) 



This is the complex conjugate of the equation for (p . The result for (pi can be obtained 
by subtracting the middle two equations in (4.13), 

= -|v/A M A r (v\ /2 (/>o + V 1/2 <h) ■ (4-17) 
The equations separate in the same way as the scalar field equation, that is we set 

<p = K (r)S ((i)e im +e- iut . (4.18) 

Then, 

(£_i/ 2 A M £j /2 - 2ux 2 M S = -XS Q 
/ / (4-19) 

(V_ 1/2 A r v\ /2 + 2iu X 2 rj n = xn . 

The second of these equations is the analogue for Kerr-de Sitter of Teukolsky's equation 
for Kerr [13]. For spin s we would have, 

(v_ s/2 A r V\ /2 + 2(2s - l)iux 2 r^ 11 = XJZ. (4.20) 

We would like to express the radial equation in terms of the r* coordinate as we 
did for the scalar field. Before doing this it is advantageous to introduce electromagnetic 
potentials rather than the fields themselves because these are directly comparable to the 
scalar field, particularly when we come to evaluate energy fluxes. In order to proceed we 
need the Teukolsky and Starobinsky identities, 

v_ 1/2 A r v l/2 n 2 = cn , vl 1/2 A^v\ /2 n = cn 2 , (4.21) 

where 

C 2 = A 2 +4u(am- a 2 u) X 4 - (4.22) 

There are similar identities for C. The derivation of these identities follows the correspond- 
ing results in reference [1]. 
Define 

R = VAl /2 n 2 , R 2 = V ] Al /2 TZ . (4.23) 
It follows from the identities that 

TZ = C^Al^VRo, TZ 2 = C~ 1 A 1 r / 2 V^R 2 . (4.24) 

Combining the two sets of radial equations gives 

V ] A r VR = CR 2 , VA r V ] R 2 = CR . (4.25) 

These expand into a form, 



(d r A r d r + K 2 /A r - 2iux 2 r) R = CR 2 
(d r A r d r + K 2 /A r + 2iu X 2 r) R 2 = CR 



(4.26) 



10 



This coupled form of the equations seems a step backwards, but it will still be acceptable 
for our analysis. The equations can be written as a scattering problem, 



d 2 

+ V s (r))a r F = (4.27) 



*2 



where 



dr 



F=[ R + R 2 \ V= {\ U , (4 _ 28) 



and 



R0-R2J' \-U V. 

_ 9 CA r A r ( rA r \ 2iurA r 
V ± = (*- m ^^--^(-^j, ^ = - W . (4.29) 

It would be possible to separate these equation using the methods described in [1], but 
this will not be necessary. 

The new radial functions can be used to define new fields, 

<S>o = R (r)S (fi)e im(l >e- i " t . (4.30) 

These are related to the electromagnetic potentials. The result (4.24) shows that the 
Maxwell fields are given by 

% = ^V^ , 02 = ^2^2, (4.31) 



and from (4.17), 



Cc ^ Cc 



0i = -\^(0o + 0i) - i^f ( $ o + $2). (4.32) 



We have not fully explored the numerous relationships that exist between the mode func- 
tions because they are not needed for our stability analysis, but there are many equivalent 
ways to express the above results. 

5. The gravitational wave equations 

In the Newman-Penrose formalism, the Weyl tensor is described by five complex 
scalars ^0 • • • ^4- Of these, only ^2 is non-vanishing in black hole solutions. The remain- 
ing scalars can be used to describe gravitational wave perturbations of the metric. In this 
section we will separate the equations for these perturbations on the Kerr-de Sitter back- 
ground. The perturbed metric has other degrees of freedom, but we leave those equations 
for another occasion. 

The Bianchi identities for the Weyl tensor components read 

(6 - Au + 7t)# - (D - 2e - 4p)tf 1 = 3k^ 2 
(A-4 7 + ^)^ -(5-4r-2/3)^i =3a^ 2 
(5 + 4(3- r)* 4 - (A + 2 7 + 4^3 = -3z/^ 2 
(D + 4e - p)* 4 - (6 + 4tt + 2a)^ 3 = -3A* 2 - 
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The connection components k, a, A and v all vanish for the unperturbed metric and the 
values appearing here are therefore the perturbed values. We need two Ricci identities to 
solve for these, 



(D-p-p-Se + e)a-(8-T + 7T-a- 3(3)k = ^ 
(A - /I - p, + 37 - 7) A - (6 + 3a + (3 + n - f )u = -* 4 . 



(5.2) 



The angular and radial operators C and V that where used for the Maxwell fields 
can be used here also. We rescale the Weyl scalars first, ipi = c 2 ^i and the connection 
components k = c x I 2 {c)~ x I 2 k. Then, 



{£{ + -) Vo 



3?: 



V[ - -_ j Vo 



3i 
c 



A M ( d + - j V4 + V&r [ Vl 1/2 + 3)^8 



3z 



£_ 1/2 - - ) ^3 



(z>_i/ 2 + Vi = 3v / 2(# 2 )/i, 
-3V2(c0 2 )(t, 

-3a/2(c0 2 )i>, 
A7 ^X>i - ^4 = 3>/2(c02)A. 



(5.3) 



and 



The Weyl scalar 



.: 3i 
£-1/2 • 7 ' 



*5 



+ ^A~ r ( Pl 1/2 + - ) A = V2c(c)- 2 if>4 



(5.4) 



M Q 2 c 



(C)^ (C)^- (5 ' 5) 

If the charge Q = 0, we can eliminate ifii from the first two equations by using equation 
(5.4). We will restrict our attention to this case alone, then 



r 3z 

L--1/2 - 



Si 



(c) 2 



After simplification, this gives 

(C-tApCl + V- X A r V\ + 6iux 2 c - 2Ap 2 ) Vo = 0. 



(5.6) 



(5.7) 



Similarly, the Weyl scalar 1^4 satisfies the complex conjugate equation. The remaining field 
ip3 is gauge dependent and can be made to satisfy similar equations or made to vanish, 
depending on the choice of tetrad in the perturbed metric. 
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We separate the equation as before with a radial function TZo(r) and an angular 
function <S (/-*). These functions satisfy 

f£_iA M £{ - Qu X 2 V - 2Afi 2 ) S = -XS 

/ / (5-8) 

\V-iA r T>{ + Qiux 2 r - 2Ar 2 J ft = Aft . 

These equations are consistent with the new version of Teukolsky's equation which was 
proposed in the previous section, but with the addition of a cosmological constant term. 
This term plays a similar role to a mass in the case of a massive scalar. 

As before, we will introduce potentials to facilitate the stability analysis. This makes 
use of the identities, 

V-tAVVAVtTli = CUo, V^AV^V^ AV\TZ = CU 4 . (5.9) 



The modulus of the constant C is given by 

| C\ 2 =\ 2 \\ - 20K K 2 \\i + 24iffa 2 A + K 2 (K^ + M 2 ) 
+ 8A(a 2 AAi - 4K 2 X + 6a 2 K K 2 ) + 16A 2 a 4 



(5.10) 



where 

Define 
and also 



K = (am-a 2 uj)x 2 , K 2 = -2u X 2 , Ai = A + 2 - |Aa 2 . (5.11) 
#0 = VVAV{R^ # 4 = V^V^AV\n , (5.12) 



$i = R l {r)S l {9)e imct, e- lujt . (5.13) 



The equations for these radial functions form a coupled set of second order equations. 
(This can be shown by repeated use of the original equation to reduce the order, but the 
resulting equations are quite complicated). The Weyl curvature scalars are recoverable 
from, 

^4 = ^2 ^^4, ^0 = ^ V $ 0- (5-14) 

The remaining scalars and ^3 can be chosen to be zero by adopting a suitable gauge. 
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6. Regularity at the Cauchy horizon 

We shall now examine the behaviour of an incoming flux of energy as it reaches the 
Cauchy horizon. The energy flux will be set up initially in the form of waves that enter 
the black hole from the outside. These waves are affected by the geometry of spacetime 
in two ways, their wavelengths are stretched or compressed and they get scattered. It is 
the compression of the waves that leads to instability, but it is necessary to consider the 
scattering problem in some detail. We begin therefore with some general remarks about 
scattering problems. The following treatment is based very closely on Chandrasekhar and 
Hartle [16]. 

Consider the following scattering problem 

^ + V(x)f = (6.1) 

where the potential has an exponential expansion as x — > ±oo, 

f + -2riK + x 

V(x) = ^° C :_ w , * °° • (6.2) 



Suppose also that Cq is parameterised by a frequency u>, 

4 = (w - mO±) 2 . (6.3) 



A potential of this form is called Yukawian [17] and leads to special analyticity properties 
of the scattering coefficients which are derived in the appendix. 

In the case of waves entering the black hole there are two scattering regions, the first 
outside the hole in region I and the second inside region II. The radial modes will be 
separated into <— modes which enter from the cosmological horizon and — > modes that are 
scattered from the past event horizon. These are shown in figure 5. Near to the Cauchy 
horizon, the asymptotic limits of the modes are 

a r R ~ Ae lkar * + Be~ lk ^ (6.4) 



in both cases. 

The scattering amplitudes can be found by combining two scattering problems with 
coordinates x = r* in both regions. This gives a potential of the required form. The 
scattering problems are related to those in the appendix by the substitutions 



I : r_ = r 2 , r + = n, f 2 = a r R, fi = a r ^; 
II : r_ = r 2 , r+ = r 3 , f 2 = (T II /T I )a r < R, h = (T n / 'R^H. 

Consequently, 



(6.5) 



(6.6) 
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The full expansion of the field $ will be of the form $ = + ~$> , with 

t = \"[ —W(lmu)*R(u;r)S(lmu;6)e im ^- tujt . (6.7) 

, J 27T 

L,m 

The functions W(lmu) determine the initial conditions on an initial surface in region I. 
We will take W to be a regular function of uj with isolated poles. 

Close to the Cauchy horizon, the asymptotic form of the mode functions can be used, 
and then 

/A 
——W{lmuj)S{lmu; 6)(Xe- ik3V + *Be ik * u ) e im ^~ Qt \ (6.8) 
27ia r 

Recall that <p — fl(rs)t is the regular coordinate fa. Futhermore, ks is taken to be real or in 
the lower half plane and therefore the exponential terms converge at the Cauchy horizon, 
where v — > oo. In fact, the only possible divergences in the field $ at the horizon would 
arise from poles in A(u), and these can be removed by taking the principal value of the 
integral. 

The energy flux in the fields would be measured by the components of the stress- 
energy tensor. A coordinate frame that is regular on the Cauchy horizon, such as the 
Kruskal coordinate frame, should be used. A typical component of the stress-energy for a 
scalar field would be 

T vv = -d v $d v $. (6.9) 

The part of the stress-energy tensor that is able to diverge at the Cauchy horizon is the 
derivative with respect to V. This is because 

d v <S> = e K3V d v $. (6.10) 

At the Cauchy horizon, v — > oo, and the exponential diverges. 

For the electromagnetic field, the stress tensor components in the null tetrad frame 
are functions of the Maxwell fields <f>o, 4>i and fa, for example 

T ab l a l b = fafa, T ab n a n b = fafa, T ab l a n b = fafa. (6.11) 

This frame has to be rescaled to obtain the Kruskal frame, 

d v = e ^("+«)/2 / 5 du = e -K 3 (u+v)/2 n _ ( 6 _ 12 ) 

The potentially diverging stress-energy tensor component is then 

T w = ^(dv*o)*(dv*o), (6.13) 
after using expresions for the Maxwell fields given in equations (4.31). 
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In the case of the gravitational waves we are interested in whether components of the 
curvature tensor remain finite at the Cauchy horizon. We have 

C abcd l a m b l c m d = -* , C abcd n a m b n c m d = -* 4 - (6.14) 

Again, the important terms are the ones depending upon the V derivatives of (The 
basis vectors m and m have to be Lorentz boosted to a regular frame.) 

Using the mode decomposition above and the asymptotic forms of the mode functions 
we can express the potentially divergent term as 




(6.15) 



and a similar term with the — > amplitude. The integrand is regular with isolated poles 
which allows us to move the contour of integration away from the real axis. The behaviour 
of the whole term as v — > oo depends upon the poles in the integrand in the lower half of 
the complex u> plane. 

The first such pole in k%A can be found from equation (8.11) and figure 6. The crucial 
feature is that the pole at u> = mVt2 — in * s matched by an identical pole in 1/Tj 

or RijTj. The ratio (6.6) therefore has no pole there. The pole at oj = m0 3 also cancels 
with k$, leaving the first pole at u> = mH 3 -jK 3 . 

In order to discuss the behavior of W we can consider in a similar way the stress- 
energy tensor component Tyy at the cosmological horizon. Here Ks is replaced by k±, 
and we see that the value of Tyy on the horizon is determined by the residue of W at 
lo = mVti — iki. Therefore, if the stress energy tensor is non-zero at the cosmological 
horizon, there should be a pole in W with imaginary part — ik\. 

We have stated already our assumption that the angular functions are analytic in a 
strip bordering the real axis. We conclude that it is possible to distort the contour of 
integration to pass parallel to the real axis through — ik\ but no further, and that Tyy 
at the Cauchy horizon diverges as exp(ft 3 — K\)v. In particular, the stress-energy tensor 
of the waves is regular at the Cauchy horizon if the surface gravity there is less than the 
surface gravity at the cosmological horizon. 

7. Conclusion 

We have seen that there are situations in which the energy fluxes of scalar, electro- 
magnetic and gravity waves remain finite at the Cauchy horizon and allow an observer to 
pass safely through into the interior of the black hole. The angular velocity of the hole 
has to be finely tuned to allow this to happen, and of course we assumed the existence of 
a cosmological constant. 

The separation of the field equations that enabled us to make progress was closely as- 
sociated with the algebraic character of the curvature tensor of the black hole and with the 
symmetrical appearance of the metric with respect to the angular and radial coordinates. 
This also explains why the angular and radial equations have such similar properties [1]. 
We have been able to predict the form of the radial field equations with a cosmological 
constant for arbitrary spin, thus generalising Teukolsky's results [13]. We have also found 
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it necessary to assume some properties of the angular mode functions that merit further 
study. 

The gravitational wave equations form a major part of the full stability analysis for this 
spacetime. We have left the complete separation of the gravitational perturbations in the 
metric for another time. The remaining metric components are mostly gauge parameters 
of little physical interest in the absence of other fields. The more interesting case would 
include electromagnetic fields as well, but separation is then problematic, even in the 
standard analysis, particularly when the hole carries an electric charge. 

Even in the case that the cosmological constant is very small, there is still a tiny 
range of values for the rotation rate of the hole leading to a stable Cauchy horizon. This 
range of values is very close to the extreme limit in which the surface gravity of the event 
horizon vanishes, a situation that we would expect only to be attainable asymptotically. 
It therefore seems possible that attempts to spin up a black hole could be used in principle 
to stabilise the Cauchy horizon. 

The stable situation raises the question of how to describe a naked singularity. The 
classical theory needs to be modified to do this, but we might make some progress by 
treating the interior geometry as a central force problem with an undetermined set of 
phase shifts at the singularity. The next question would be what are the consequences of 
closed timelike curves in the interior of the hole. There is a possibility that linear wave 
propagation and therefore the stability analysis could still be consistently defined, as they 
are in asymtotically flat spacetimes [18]. Finally, we would be left with the continuation 
of the observer's journey through the hole and into another universe. 

8. Appendix 

In this appendix we describe the procedure by which the poles of the scattering am- 
plitudes for equation (6.1) can be found, following ref [16] but allowing the potential to 
approach different constant values at the horizons. Define solutions of the scattering prob- 
lem fi and f 2 with the asymptotic behaviour 

/ 1 (x, W )-e +ifc + !B h{x,u>) ^ e~*+* x^oc 

f 2 (x, uj) -> e ~ ik - x f 2 (x, uj) -> e +ik ~ x x -> -oo 

where 

k±=u)-mn±. (8.2) 

(Note that ± refer to ±oo, the reverse of ref. [16].) 

The solutions fi and fi form a complete set of solutions to the scattering problem 
and it is therefore possible to express f 2 as a linear combination, 

f 2 (x,uj) = ^If^x.oj) + -L-ftfau). (8.3) 

1 [UJ) 1 [UJ) 

This defines the reflection and transmission coefficients R(u) and T(oj). They can be 
extracted from fi and f 2 for any value of x by using the Wronskians, 

■[fi(x,uj),f 2 (x,u)] 



Tiuo) 2ik+ 

R(J) 1 ' (8 - 4) 

^^- [/l (, W ),/ 2 (, W )] 
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where [/, g] = fg' - gf. 

The scattering problem can be reduced to an integral equation by using the Green 
function. We first write the differential equation as 

2 • /•-'/,' OV) r )/..,. (8.5) 



dx 2 

The solution obtained by regarding the right hand side as a source term is 

r x 1 

f 2 (x,u) = e~ lk - x + / —smk-{x-x')(V{x')-CQ)f 2 {x',u>)dx'. (8.6) 
The integral equation can be iterated to obtain a series for f 2 . set 

oo 

f 2 (x,u) = e-* k -* + Y,f2 S) (x^) (8.7) 

s=l 

then, 



/x y 
— sin/c_(x - x s ) (V(x s ) - eg) f 2 (x s , u) dx s . 

Repeatedly iterating this equation gives 



(8.8) 



oo 

2 ?7' At — X n 

/■■ 

n 



(8.9) 

where the exponential expansion of V has been inserted. If we define j/j = — a^, then 
following [16], we can write 



,( S ), e 



-ik-x r°° 



n s, (x 1 oj) = — 



[2ik-Y 



/ dy\ • • • I dy s H\(l-e 2ik -yi)J2cne 2nK - Xi \- (8-10) 
J ° Jo i=1 [ n=1 J 



The integrand is a sum of exponentials and the result of integrating over the yi intro- 
duces poles at k- = —in,K-, where n = 1,2, . . .. Therefore, f 2 (x, ui) has poles at 

u> = mfi_ — inK-. (8-11) 

A similar treatment of fi(x,u) shows that this function has its poles at 

u> = mO + + inK + . (8-12) 

From equation (8.4), these determine where the poles of R(uj)/T(u) and 1/T(u>) lie, as 
shown in figure 6. 
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Figures 

1. The parameter space of the black hole for (a) Q = and (b) Q = 0.2. The mass M 
and rotation rate a have both been scaled by a, where a 2 = 3/A. Coincident horizons 
occur along the outer borders of the triangular region, t2 = rs on OA and r± = r2 
on AC. Corresponding surface gravities also vanish there. Interior lines denote the 
cases where k± = Ks (the lower line OA) and k± = K2 (along OB). The narrow region 
extending from O to A is the stability region. 

2. A Penrose diagram showing the event horizon T2 between regions I (outside the hole) 
and II (inside the hole). 

3. Part of the Penrose diagram of the fully extended spacetime. Horizon 1 is a cosmo- 
logical de Sitter horizon, horizon 2 a black hole event horizon and horizon 3 an inner 
Cauchy horizon. Asymptotic regions exist and are labelled as J. Singularities at r = 
are represented by zigzag lines. The fully extended diagram covers the whole plane, 
including the (blank) regions that can only be reached by passing through r = 0. 

4. The radial potential (reversed sign) in region II for various values of /. The upper 
three curves have m = 1 and the lower three m = 0. 

5. The inward and outward radial mode functions R and ^ showing how their asymp- 
totic amplitudes A, B and C are associated with the respective horizons. 

6. The location of the poles in 1/T and R/T for s = are shown in the complex w-plane. 
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Tables 



lo a 


(Z,m) = (0,0) 


(1,0) 


(1,1) 


(1,-1) 


(2,0) 


0.0 


0.000 


2.000 


2.000 


2.000 


6.000 


0.1 


0.007 


2.004 


1.809 


2.208 


6.005 


u.z 


0.027 


2.016 


1.631 




« ni q 


3 


0.060 


2.036 


1.471 


2 672 


6 043 


4 


0.107 


2.064 


1.327 


2.982 


6 076 


1 


0.667 


2.400 


0.800 


4 800 


6 476 




Table 1. Angular eigenvalues for a/ 


a = 0. 




u> a 


(Z,m) = (0,0) 


(1,0) 


(1,1) 


(1,-1) 


(2,0) 


0.0 


0.000 


2.003 


2.021 


2.097 


6.208 


0.1 


0.007 


2.008 


1.825 


2.310 


6.212 


0.2 


0.028 


2.020 


1.645 


2.538 


6.228 


0.3 


0.061 


2.040 


1.482 


2.783 


6.252 


0.4 


0.109 


2.069 


1.335 


3.044 


6.249 


1.0 


0.680 


2.473 


0.857 


4.954 


6.694 



Table 2. Angular eigenvalues for a/ a = 0.1. 



u; a 


(Z,m) = (0,0) 


(1,0) 


(1,1) 


(1,-1) 


(2,0) 


0.0 


0.000 


2.007 


2.079 


2.403 


6.861 


0.1 


0.008 


2.012 


1.871 


2.628 


6.867 


0.2 


0.029 


2.025 


1.670 


2.871 


6.881 


0.3 


0.065 


2.046 


1.508 


3.130 


6.906 


0.4 


0.116 


2.077 


1.352 


3.407 


6.943 


1.0 


0.721 


2.700 


1.040 


5.432 


7.375 



Table 3. Angular eigenvalues for a/a = 0.2. 
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